Generation and Propagation of Nonlinear Waves in Travelling Wave Tubes 
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The generation and evolution of nonlinear waves in microwave amplifiers such as travelling wave 
tubes, free electron lasers and klystrons have been studied. The analysis is based on the hydro- 
dynamic and field equations for the self-consistent evolution of the beam density distribution, the 
current velocity and the electromagnetic fields. A system of coupled nonlinear Schrodinger equa- 
tions for the slowly varying amplitudes of interacting beam-density waves has been derived. Under 
the approximation of an isolated mode neglecting the effect of the rest of the modes, this system 
reduces to a single nonlinear Schrodinger equation for that particular mode. 
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I. INTRODUCTION 

The generation and evolution of nonlinear waves and 
harmonic distortions in microwave amplifiers such as 
travelling wave tubes, free electron lasers (FELs) and 
klystrons have recently attracted much research interest 
In connection with the construction and com- 
missioning of the next generation of FELs and powerful 
klystrons for accelerating RF cavities in circular machines 
and linear colliders, this issue has become even more chal- 
lenging. Of particular importance are the effects of in- 
tense self- fields due to space charge and current, as well 
as wake fields due to interaction impedances. Both of 
the above influence the propagation of the electron beam 
in microwave devices, its stability and transport prop- 
erties. In general, a complete description of collective 
processes in intense charged particle beams is provided 
by the Vlasov-Maxwell equations for the self-consistent 
evolution of the beam distribution function and the elec- 
tromagnetic fields. Usually, the electron beam in a trav- 
elling wave tube can be assumed to be weakly colli- 
sional. Hence, the dynamics of electrons is well de- 
scribed by the hydrodynamic equations coupled with the 
equations for the electromagnetic self-fields, which con- 
stitutes a substantial simplification of the model. Al- 
though the analytical basis for modelling the dynamics 
and behaviour of space-charge-dominated beams is well 
established, a thorough and satisfactory understanding 
of collective processes, detailed equilibrium and forma- 
tion of patterns and coherent structures is far from being 
complete. 

While the linear theory of wave generation in mi- 
crowave amplifiers is relatively well understood £j, the 
nonlinear regime is far from being exhaustively studied. 
The present paper is aimed at filling this gap. We will 
be mainly interested in describing the slow evolution of 
some coarse-grained quantities that are easily measur- 



able, such as the wave amplitudes. Owing to the nonlin- 
ear wave interaction contingent on the nonlinear coupling 
between the hydrodynamic and Maxwell equations, one 
can expect a formation of nontrivial coherent structure 
that might be fairly stable in space and time Here, 
we show that solitary wave patterns in the electron beam 
density distribution are an irrevocable feature, character- 
istic of powerful microwave amplifiers. 

The paper is organized as follows. In the next section, 
we state the basic equations which will be the subject 
of the renormalization group (RG) reduction in section 
III. Starting from a single equation [see equation 120(1 ] 
for the density distribution of the electron beam, we ob- 
tain a formal perturbation expansion of its solution to 
second order. As expected, it contains secular terms pro- 
portional to powers of the time variable which is the only 
renormalization parameter adopted in our approach. In 
section IV, the arbitrary constant amplitudes of the per- 
turbation expansion are renormalized such as to elimi- 
nate the secular terms. As a result, a set of equations for 
the renormalized slowly varying amplitudes is obtained, 
known as the renormalization group equations (RGEs). 
These equations comprise an infinite system of coupled 
nonlinear Schrodinger equations. Finally, section V is 
dedicated to discussion and conclusions. 



II. FORMULATION OF THE PROBLEM AND 
BASIC EQUATIONS 

The electron beam in a travelling wave tube is assumed 
to be weakly collisional. Therefore, the dynamics of elec- 
trons is well described by the hydrodynamic equations 
coupled with the equations for the electromagnetic self- 
fields. We start with the ID description of a beam of 
electrons propagating in an external focusing electric field 
with focusing coefficient 67, which models the bunching 
of the electron beam in the longitudinal direction. As we 
will see in the sequel its additional role is to attain a sta- 
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tionary equilibrium by partially compensating the space- 
charge defocusing. The continuity and the momentum 
balance equations can be written as 



9q 9 , . 
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where g and v are the electron density and the current 
velocity. Furthermore, m, e and T are the mass, the 
charge and the temperature, respectively, while ks is the 
Boltzmann constant. The space-charge tp sc and wave ip w 
potentials satisfy the Poisson and the wave equation 
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where □ denotes the well-known d'Alembert operator. 
In addition, Eq is the permittivity of free space and 
n = N e /Vt is the electron number density, where N e 
is the total number of electrons and V t is the volume 
occupied by the electron beam in the longitudinal direc- 
tion. Moreover, the electron beam cross-sectional area is 
represented by A, the quantity Zq denotes the interac- 
tion impedance, while c represents the phase velocity of 
a cold circuit wave. 

Let us introduce the scaling 



en 



¥>W = W, 



and rewrite equations as follows 
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equation JSJ can be alternatively expressed according to 
the relation 



v T = uj p r D , 



(10) 



where td is the electron Debye radius. Equations jTfr. 
© and J7J) possess a stationary equilibrium solution 
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Therefore, we can further scale the hydrodynamic and 
field variables as 

g = g Q + eR, U sc = U + ell, v -> ev, W -> eW, 

(12) 

where e is a formal small parameter introduced for con- 
venience, which will be set equal to one at the end of the 
calculations. Thus, the basic equations to be used for the 
subsequent analysis can be written in the form 
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Before we continue with the renormalization group re- 
duction of the system of equations 1131) 1 51) in the next 
section, let us assume that the actual dependence of the 
quantities i?, v, U and W on the spatial variables is rep- 
resented by the expression 



* = V = (R,v,U,W), 



(16) 



where £ = ez is a slow spatial variable. Thus, the only 
renormalization parameter left at our disposal is the time 
t which will prove extremely convenient and simplify te- 
dious algebra in the sequel. 



The electron plasma frequency u> p and the thermal ve- 
locity vt of the electron beam are expressed according to 



III. RENORMALIZATION GROUP 
REDUCTION OF THE HYDRODYNAMIC 
EQUATIONS 
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(8) Following the standard procedure of the renormaliza- 

tion group method, we represent f as a perturbation 
expansion 
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is a shorthand parameter introduced for later conve- 
nience. Note that the thermal velocity vt as defined by 



in the formal small parameter e. The next step consists 
in expanding the system of hydrodynamic and field equa- 
tions l|13|) - (|15(l in the small parameter e, and obtaining 
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their naive perturbation solution order by order. It is 
possible to simplify this system, which will turn out ex- 
tremely useful in what follows. Differentiating equation 
(|13[l with respect to the time t, differentiating equation 
(|14|l with respect to z and using equations JTSJ, we can 
eliminate the electric potentials. As a result of obvious 
manipulations, we obtain 
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where the linear differential operator C(d z ,dt) is given 
by the expression 
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Taking into account the expansion of the logarithm, we 
notice that the right-hand-side of equation (|18fl is at least 
of first order in the formal parameter, so that 
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Equation (|20|l represents the starting point for the renor- 
malization group reduction, the final goal of which is to 
obtain a description of the relatively slow dynamics lead- 
ing to formation of patterns and coherent structures. 

Let us proceed order by order. The solution to the 
zero-order perturbation equation i|20|) can be written as 



iZo(*,£;t)=X>k(0< 
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and A k is an infinite set of constant complex amplitudes, 
which will be the subject of the renormalization proce- 
dure in the sequel. Here "constant" means that the am- 
plitudes Ak do not depend on the fast spatial variable 
z and on the time t, however, they can depend on the 
slow spatial variables £. The summation sign in equation 
(|21|l and throughout the paper implies summation over 
the wave number k in the case where it takes discrete 
values, or integration in the continuous case. From the 
dispersion equation 
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it follows that the wave frequency to k can be expressed 
in terms of the wave number fc, where the Fourier-image 



□fe of the d'Alembert operator can be written according 
to 



It is important to emphasize that 
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where the asterisk denotes complex conjugation. The 
latter assures that the perturbed density distribution as 
defined by equation (|21|l is a real quantity. The zero- 
order current velocity vo(z 1 S l ;t) obtained directly from 
equation H13|l can be written as 
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In first order equation (|20|l acquires the form 
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where by C z we have denoted the derivative of the oper- 
ator C with respect to d z . It has now two types of solu- 
tions. The first is a secular solution linearly dependent on 
the time variable in the first-order approximation. As a 
rule, the highest power in the renormalization parameter 
of the secular terms contained in the standard pertur- 
bation expansion is equal to the corresponding order in 
the small perturbation parameter. The second solution 
of equation (|20|l arising from the nonlinear interaction 
between waves in the first order, is regular. Taking into 
account the fact that the Fourier image of the operator 
C z is equal to —i(dT>/dk), we can write the equation for 
determining of the secular part of the solution as 
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dujk <Juj k 
which is a direct consequence of the general relation 

CtG(t) = tCG(t) + C t G(t), (30) 

holding for a generic function G(t). Here C t implies dif- 
ferentiation with respect to <9 t . To verify equation l|30|l . 
it suffices to prove the identity by induction in the case, 
where C is the monomial operator 9™ and then take 
into account the Taylor expansion of C. With these re- 
marks in hand, it is straightforward to solve equation 
(J2BJ. Combining its solution with the solution of the 
regular part, we obtain 
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where u g k is the group velocity defined as 
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In explicit form, the components of the infinite matrix 
ctki are given by the expression 
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Furthermore, the first-order current velocity can be ex- 
pressed as 
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where 



second-order perturbation equation (|20|l will contribute 
to the secular solution. Since we are interested in this 
solution to be renormalized later, the contributing terms 
will be retained only. Hence, we can write 



r p + ? dRl 4. c " d2Ro n 



£>o- 



9 2 _(woWi) 
<9z 2 



dtdz 

or in explicit form 

^2 = E 



(40) 



dV d 2 V \ ld 2 V 



dk dkdu k J 2 dk 2 



d 2 A k 



'Jo 



(41) 



k,m 



The components of the infinite matrix T mk can be ex- 
pressed according to the relation 
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A standard algebra similar to the one outlined in the 
previous section leads to the second order secular solution 
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A couple of interesting features of the zero and first- 
order perturbation solution are noteworthy to be com- 
mented at this point. First of all, the zero-order den- 
sity, current velocity (and electric potentials) are pro- 
portional to the arbitrary complex amplitudes A k . The 
second terms in the expressions for the first-order den- 
sity Ri and current velocity v\ [see equations (|31f) and 
(|37|l ] imply contribution from nonlinear interaction be- 
tween waves. It will be shown in the remainder that these 
terms give rise to nonlinear terms in the renormalization 
group equation and describe solitary wave behaviour of 
a generic mode. 



IV. THE RENORMALIZATION GROUP 
EQUATION 

Passing over to the final stage of our renormalization 
jroup procedure, we note that particular terms in the 
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Taking into account the definition l|32|) of the group ve- 
locity, we conclude that 
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Following the standard procedure [U |(| of the RG 
method, we finally obtain the desired RG equation 
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where now A k is the renormalized complex amplitude p| . 
Thus, the renormalized solution for the density pertur- 
bation of the electron beam acquires the form 

R(z-t) = Y J Mz;t)^ k(z '' t \ (47) 

k 

For the renormalized electric field, we obtain 

E(z; t) = £ fc24 ~^ fc (z; t)e^\ (48) 
e G ^ fc 

Equations (|46[) represent a system of coupled nonlin- 
ear Schrodinger equations for the amplitudes of eigen- 
modes. Consider a particular mode with a wave number 
k. As a first approximation, the contribution of the other 
modes with m ^ k can be neglected, which results in a 
single nonlinear Schrodinger equation for mode k. The 
nonlinearity in the corresponding nonlinear Schrodinger 
equation describes the nonlinear interaction of the mode 
k with itself. 

It should be emphasized that the approach outlined 
in the present paper is rather general even in the the 
case of more than one renormalization parameter, where 
the extension should be straightforward. Note that the 
unique assumption concerns the knowledge of the dis- 
persion properties of the linear differential operator gov- 
erning the evolution of the zero-order quantities. The 



approach can be successfully applied to a wide class of 
problems of particular physical interest. 

V. DISCUSSION AND CONCLUSIONS 

We have studied the generation and evolution of non- 
linear waves in microwave amplifiers such as travelling 
wave tubes, free electron lasers and klystrons. The anal- 
ysis performed in the present paper is based on the hydro- 
dynamic and field equations for the self-consistent evo- 
lution of the beam density distribution, the current ve- 
locity and the electromagnetic fields. Using further the 
RG method, a system of coupled nonlinear Schrodinger 
equations for the slowly varying amplitudes of interact- 
ing beam-density waves has been derived. Under the 
approximation of an isolated mode neglecting the effect 
of the rest of the modes, this system reduces to a single 
nonlinear Schrodinger equation for that particular mode. 

Since the approach pursued here is rather general, it is 
presumed that it may find applications to other problems, 
where of particular interest is the dynamics of slowly 
varying amplitudes of patterns and coherent structures. 
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